Exact results for curvature-driven coarsening in two dimensions 



Jeferson J. Arenzon, 1 Alan J. Bray, 2 Leticia F. Cugliandolo, 3 and Alberto Sicilia 3 

'instituto de Fisica, Universidade Federal do Rio Grande do Sul, CP 15051, 91501-970 Porto Alegre RS, Brazil 

^School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK 
Universite Pierre et Marie Curie - Paris VI, LPTHE UMR 7589, 4 Place Jussieu, 75252 Paris Cedex 05, France 

(Dated: February 6, 2008) 

We consider the statistics of the areas enclosed by domain boundaries ('hulls') during the 
curvature-driven coarsening dynamics of a two-dimensional nonconserved scalar field from a dis- 
ordered initial state. We show that the number of hulls per unit area that enclose an area greater 
than A has, for large time t, the scaling form Nh(A,t) = 2c/ (A + At), demonstrating the validity of 
dynamical scaling in this system, where c = l/87T\/3 is a universal constant. Domain areas (regions 
of aligned spins) have a similar distribution up to very large values of A/Xt. Identical forms are 
obtained for coarsening from a critical initial state, but with c replaced by c/2. 



Coarsening dynamics has attracted enormous interest 
over the last 40 years. The classic scenario concerns a 
system that in equilibrium exhibits a phase transition 
from a disordered high-temperature phase to an ordered 
low-temperature phase with a broken symmetry of the 
high-temperature phase. The simplest example is, per- 
haps, the Ising ferromagnet. When the system is cooled 
rapidly through the transition temperature, domains of 
the two ordered phases form and grow ('coarsen') with 
time under the influence of the interfacial surface ten- 
sion, which acts as a driving force for the domain growth 

ail. 

While phase transitions provide the traditional arena 
for coarsening dynamics, there are many other examples, 
including soap froths Q , breath figures [f| , granular me- 
dia @, and interfacial fluctuations Q. A common fea- 
ture of nearly all such coarsening systems is that they are 
well described by a dynamical scaling phenomenology in 
which there is a single characteristic length scale, R(t), 
which grows with time. If dynamical scaling holds, the 
domain morphology is statistically the same at all times 
when all lengths arc measured in units of R(t). The as- 
sumption of dynamical scaling also makes possible the 
determination of the length scale R(t) for a large class of 
coarsening systems 0, U 

Despite the success of the scaling hypothesis in de- 
scribing experimental and simulation data, its validity 
has only been proved for very simple models, includ- 
ing the Id Glauber-Ising model 0] and the nonconserved 
0(n) model in the limit n — > oo [Iol |. Another note- 
worthy exact result is the Lifshitz-Slyozov derivation of 
the domain-size distribution for a conserved scalar field in 
the limit where the minority phase occupies a vanishingly 
small volume fraction The only other exact results, 
to our knowledge, for domain-size distributions in coars- 
ening dynamics are for the zero-temperature Glauber- 
Potts and time-dependent Ginzburg-Landau [HI 
models in Id. 

In the present work we obtain some exact results for 
the coarsening dynamics of a nonconserved scalar field in 
d = 2, demonstrating explicitly, en passant, the validity 



of the scaling hypothesis. To do this, we use a contin- 
uum model in which the velocity, v, of each element of a 
domain boundary is proportional to the local interfacial 
curvature, k: 

v = -(A/2tt)k , (1) 

where A is a material constant with the dimensions of a 
diffusion constant, and the factor 1/2tt is for later conve- 
nience. The Allen-Cahn equation (|TJ) may be derived 
from the zero-temperature time-dependent Ginzburg- 
Landau equation for the underlying order-parameter field 

i!. 

From Eq. |T]), we can immediately deduce the time- 
dependence of the area contained within any finite hull 
(i.e., the interior of a domain boundary) by integrat- 
ing the velocity around the hull: dA/dt = §vdl — 
— (A/27r) § ndl = —A, the final equality following from 
the Gauss-Bonnet theorem. At any given time t, there- 
fore, hulls with original enclosed-area smaller than At 
will have disappeared, and the enclosed-areas of surviv- 
ing hulls will have decreased by Xt. In other words, the 
entire distribution of hull enclosed areas is advected uni- 
formly to the left at rate A. If Nh(A,t) s the number 
of hulls per unit area of the system with enclosed-area 
greater than A, it follows that 

N h {A,t) = N h {A + Xt,0) , VAX). (2) 

To determine the initial condition we note that, shortly 
after the quench from the high temperature phase, the 
system is at the critical point of continuum percolation. 
Cardy and Ziff [3| have shown that the number of perco- 
lation hulls, per unit area of the system, with area greater 
than A has, for large A, the universal asymptotic form 

N p (A) ~ c/A (3) 

where c = l/87r-\/3 is a universal constant. This re- 
sult provides the desired initial condition, Nh(A,0) — 
2N p (A), in Eq. ©, giving 

N h (A,t) = 2c/(A + Xt) , (4) 



2 



where the factor 2 arises from fact that there are two 
types of hull, corresponding to the two phases, while 
the Cardy-Ziff result accounts only for clusters of occu- 
pied sites (and not clusters of unoccupied sites). From 
this result one immediately derives the hull enclosed 
area density function, rih(A,t) = —dNh(A,t)/dA, where 
rift (A, t) dA is the number of hulls, per unit area of the 
system, having area in the interval (A, A + dA) : 



n h (A,t) = 2c/(A + Xt) 2 



(5) 



Equation {4} has the expected scaling form Nh(A, t) = 
t^ 1 f(A/t) corresponding to a system with characteristic 
area proportional to t. This corresponds to character- 
istic length scale R(t) ~ t?, which is the known result 
if scaling is assumed Here, however, we do not as- 
sume scaling - rather, it emerges from the calculation. 
Furthermore, the conventional scaling phenomenology is 
restricted to the 'scaling limit': A — > oo, t — > oo with 
A/t fixed. Equation by contrast, is valid whenever 
t is sufficiently large, and does not (at least on the con- 
tinuum) require large A. This follows from the fact that, 
for large t, the form Q probes, for any A, the tail (i.e. 
the large- A regime) of the Cardy-Ziff result (J3J), which is 
just the regime in which the latter is valid. 

It is, however, instructive to consider what can be de- 
duced from scaling alone, augmented by the drift equa- 
tion |2|). The general scaling form corresponding to a 
scale area t is Nh(A,t) — f(A/t), with arbitrary ex- 
ponent 4>. Consistency with requires that Nh{A,t) 
depends on A only through the combination A + At, forc- 
ing N h (A, t) cx (A + At)"^. Finally, cf> = 1 is fixed by 
the requirement that there be of order one hull per scale 
area, i.e. N h (0,t) cx 1/At. Thus, for an internally con- 
sistent picture, one requires Nh(A, 0) cx A^ 1 for large A, 
and it is gratifying that the Cardy-Ziff result not only 
has this form but also provides the exact value of the 
proportionality constant. 

The argument above relies on the T = Allen-Calm 
Eq. (JXJ) . Temperature fluctuations have a two-fold effect. 
On the one hand they generate equilibrium thermal do- 
mains that are not related to the coarsening process. On 
the other hand they roughen the domain walls thus op- 
posing the curvature driven growth and slowing it down. 
Once equilibrium thermal fluctuations are subtracted - 
equivalently, hulls associated to the coarsening process 
are correctly identified - the full temperature dependence 
should enter only through the value of the parameter A, 
which sets the time scale. For simplicity we focus here 
on zero working temperature. In a future publication we 
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shall show the finite T effects discussed above 

To test the above result we carried out numerical sim- 
ulations on the 2d square-lattice Ising model (2tflM) with 
periodic boundary conditions using a heat-bath algo- 
rithm with random sequential updates. All data have 
been obtained using systems with size L 2 = 10 3 x 10 3 
and 2 x 10 3 runs using independent initial conditions. 
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FIG. 1: (Colour online) Number density of hulls per unit 
area for the zero temperature dynamics of the 2d Ising model 
evolving from an infinite temperature initial condition. The 
full line is the prediction ((5} with c = l/87r\/3 and A = 2.1. 
Inset: finite size effects at t — 16 MCs; four linear sizes of the 
sample are used and indicated by the datapoints. The value of 
A/t at which the data separate from the master curve grows 
very fast with L with an exponent close to 2. 



We mimic an instantaneous quench from infinite temper- 
ature with a random initial condition with spins pointing 
up or down with probability 1/2. The data are plotted 
in log-log form to test the prediction nh(A, t) cx A~ 2 for 
large A. The data are in remarkably good agreement with 
the prediction ([5]) - shown as a continuous curve in Fig. 1 
- over the whole range of A and t. The downward devia- 
tions from the scaling curve are due to finite size effects. 
The latter are shown in more detail in the inset where we 
display the t — 16 MCs reults for four linear sizes. Finite 
size effects appear only when the weight of the distri- 
bution has fallen by many orders of magnitude (7 for a 
system with L = 10 3 ) and are thus quite irrelevant. The 
only fitting parameter is A, which has the value A = 2.1 in 
Fig. 1. The agreement between theory and data is all the 
more impressive given that the curvature-driven growth 
underlying the prediction (|5|) only holds in a statistical 
sense for the lattice Ising model |1 51 ] . 

The mean hull enclosed area, per unit area of 
the system, diverges with the system size: (A) = 

r o L dAAn h (A,t) ~ 2cln(i 2 /At). The fact that the total 
hull enclosed area exceeds the area of the system seems 
paradoxical until one notes that a given point in space 
belongs to many hulls. 

It is clear that the evolution of the hull-enclosed area 
distribution follows the same 'advection law' ([2]), with 
the same value of A, for other initial conditions. More- 
over, Eq. §5§ applies to any Tq > T c equilibrium initial 
condition asymptotically. Equilibrium initial conditions 
at different To > T c show only a different transient be- 
haviour; initial states that are closer to T c take longer to 
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reach the asymptotic law ([2]) (see the inset in Fig. 2) [l7j |. 

An equilibrium state at the critical temperature, Tq = 
T c , is, as expected, different. This case has already been 
addressed, for general space dimension, in the context of 
coarsening from an initial state with long-ranged spatial 
correlations [2(|. It was argued that the characteristic 
scale R(t) still grows as i 1 / 2 , since the growth is curva- 
ture driven, but the space and time correlation functions 
are modified if the spatial correlations in the initial state 
are sufficiently long-ranged [20]. The statistics of hull 
and domain areas, however, were not discussed. In place 
of the continuum percolation statistics that characterize 
the initial state shortly after a quench from high temper- 
ature, the initial state for a quench from T c is character- 
ized by the statistics of Ising cluster hulls at the critical 
point. The area distribution of the hulls of these clusters 
has been studied by Cardy and Ziff [lj] . The result has a 
form identical to Eq. ([3]), but with c replaced by c/2. We 
predict, therefore, that our results for the hull-enclosed 
area distributions can be generalized to this case by sim- 
ply making the replacement c — > c/2 everywhere, and 
keeping the value of A unchanged. In Fig. [3] we compare 
this prediction to the number density of hull-enclosed ar- 
eas in the 2dIM evolving at T = from a critical initial 
condition, once more obtaining excellent agreement. 
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FIG. 2: (Colour online) Number density of hulls per unit area 
for the zero-temperature 2dIM evolving from critical initial 
conditions. We obtained the initial states after running 10 3 
Swendsen-Wang algorithm steps. The full (red) line repre- 
sents ((5} with c — > c/2 and A = 2.1. The dashed (blue) line is 
© with c. Inset: comparison between the hull-enclosed area 
distribution at t = 128 MCs for equilibrated initial conditions 
at To = 2.5 and T -> oo. 

Also interesting is the distribution of domain ar- 
eas, rid(A, t), which are the areas of regions of aligned 
spins [1 61 ] . Domains are obtained from hulls by removing 
any interior hulls. The domain area distribution rid(A, t) 
(number density of domains with area A, per unit area) 
of the 2dIM evolving at zero temperature after a quench 
from infinite temperature is shown in Fig. 3. The down- 



ward deviations from the scaling curve in the main panel 
as well as the bumps on the tail are finite size effects due 
to the percolating clusters and, as for the hull-enclosed 
areas (see Fig. 1), are not very important. 
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FIG. 3: (Colour online) Number density of domains per unit 
area for the zero-temperature 2dIM evolving from To — > oo 
(upper panel) and To = T c (lower panel) initial conditions. In 
the main panels the percolating domains have been extracted 
from the analysis while in the insets we show the same data 
including the percolating clusters. We obtained the initial 
states after running 10 s Swendsen-Wang algorithm steps. The 
full (red) line represents ((5j) with c = l/8ny/3 (upper) and 
c — > c/2 (lower), and A = 2.1 in both cases. The dotted 
(blue) lines have slopes -2.03 (a) and -2.05 (b). 

Remarkably, the domain area distribution, rid, seems 
to be almost identical to the hull-enclosed area distribu- 
tion rih, i.e.- rid ~ 2crf/(A + Xdt) T with a prefactor 2cd, 
a parameter Xd and an exponent r taking approximately 
the same values as for the hull-enclosed areas. An ar- 
gument that treats interior domain walls in a mean-field 
approximation and uses the exact result for rih derived 
above, allows one to derive [ijj 

n d (A, t) ~ 2c d [\ d (t + t )} T - 2 /[A + \ d (t + t )Y , (6) 

with \ d = A + O(c), c d = c + 0(c 2 ) and r - 2 + 0{c) for 
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infinite temperature initial conditions and, Cd — > Qj/2 for 
initial conditions equilibrated at To = T c . to is such that 
Arfi = a 2 , with a a microscopic length-scale, and sets the 
time scale. The constant c c i and r characterize the initial 
distribution of domain areas, nd(A,0) ~ 2cda T ~ 2 /U r . 
The exponent t is known for the 2dIM critical geomet- 
rical clusters, t = 379/187 ~ 2.03 [ijj], and for cluster 
masses in 2d percolation, r = 187/91 « 2.05 [3]; the 
constant has not been computed analytically, t > 2 
allows one to satisfy the sum rule establishing that the 
total domain area, per unit area of the system, must be 
equal to unity since every spin belongs to one, and only 
one, domain. Unfortunately, it is hard to distinguish a 
power law with r = 2 or 2.03-2.05 from our numerical 
data: both exponents describe the power-law tail equally 
well, as shown in Fig. [3J with a (blue) dotted line. It is, 
however, possible to put to the numerical test the value 
of the prefactor a (or 2cd) and Xd- We shall present this 
analysis elsewhere [17] . 

Interestingly, similar results are obtained for the 2d 
random ferromagnet, at least when activation is not too 
important. In Fig. 0] we display data obtained at T = 0.4 
after a quench from To — > oo. This low but finite working 
temperature is enough to avoid the complete pinning of 
domain walls by quenched disorder. Good scaling with 
the typical domain area R 2 (t) is obtained for A/R 2 (t) > 
10 _1 and the master curve resembles strongly the one 
found in the pure Ising case also included in Fig. [4] 
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FIG. 4: (Colour online) Number density of hull-enclosed and 
domain (lower inset) areas per unit area for the 2d random 
ferromagnetic Ising model with a uniform distribution of ex- 
changes in the interval [2 — e, 2 + e] with e = 0, 0.5, 1, 1.5, 2. 
The system evolves at T = 0.4 after a quench from To — * oo. 
Upper inset: the evolution of R(e, t) used to scale the pdfs 
and its comparison to the pure R(0, t) oc t 1 / 2 law. 

In summary, we have obtained exact results for the 
statistics of hull-enclosed areas during the phase-ordering 
dynamics of 2d systems coarsening under curvature 
driven growth. Notably, these results include a proof 



of the scaling hypothesis for these systems. The results 
are strongly supported by simulations of Ising systems, 
suggesting a strong degree of universality. The domain 
area distribution satifies a very similar law. 

AJB and LFC thank R. Blumcnfcld, J. Cardy, D. 
Dhar, C. Godreche, S. Majumdar, P. Sollich, and R. 
Stinchcombe for very useful discussions. AJB and LFC 
thank the Isaac Newton Institute, Cambridge, and JJA 
the LPTHE Jussieu, Paris, for their hospitality during 
the period when this work was being carried out. JJA, 
LFC and AS acknowledge financial support from Capes- 
Cofecub research grant 448/04. JJA is partially sup- 
ported by CNPq and FAPERGS and LFC by IUF. 



[1] I. M. Lifshitz, Zh. Eksp. Teor. Fiz. 42, 1354 (1962) [Sov. 

Phys. JETP 15, 939 (1962)]. 
[2] S. M. Allen and J. W. Cahn, Acta Metall. 27, 1085 

(1979). 

[3] A. J. Bray, Adv. Phys. 43, 357 (1994), and references 
therein. 

[4] W. Y. Tarn, A. D. Rutenberg, B. P. Vollmayr-Lee, and K. 

Y. Szeto, Europhys. Lett. 51, 223 (2000), and references 

therein. G. L. Thomas, R. M. C. de Almeida, and F. 

Graner, Phys. Rev. E 74, 021407 (2006). 
[5] M. Marcos-Martin, D. Beysens, J. P. Bouchaud, C. 

Godreche, and I. Yekutieli, Physics A 214, 396 (1995). 
[6] I. Aranson and L. S. Tsimring, Rev. Mod. Phys. 78, 641 

(2006). 

[7] M. Siegert, Phys. Rev. Lett. 73, 1517 (1994); L. Gol- 
ubovic, Phys. Rev. Lett. 78, 90 (1997); D. Moldovan and 
L. Golubovic, Phys. Rev. E 61, 6190 (2000). 

[8] A. J. Bray and A. D. Rutenberg, Phys. Rev. E 49, R27 
(1994); Phys. Rev. E 51, 5499 (1995). 

[9] J. G. Amar and F. Family, Phys. Rev. A 41, 3258 (1990); 
A. J. Bray, J. Phys. A 23, L67 (1990). 
[10] A. Coniglio and M. Zannetti, Europhys. Lett. 10, 575 
(1989). 

[11] I. M. Lifshitz and V. V. Slyozov, J. Phys. Chem. Solids 
19, 35 (1961). 

[12] B. Derrida and R. Zeitak, Phys. Rev. E 54, 2513 (1996). 
[13] A. J. Bray, B. Derrida, and C. Godreche, Europhys. Lett. 

27, 175 (1994). 
[14] J. Cardy and R. M. Ziff, J. Stat. Phys. 110, 1 (2003). 
[15] A. D. Rutenberg, Phys. Rev. E 54, R2181 (1996) has 

shown that lattice anisotropy also induces a very small 

anisotropy in the coarsening structure. 
[16] See T. B. Liverpool, Physica A 224, 589 (1996) and 

M. Fialkowski and R. Holyst, Phys. Rev. E 66, 046121 

(2002) for earlier studies of this problem. 
[17] J. J. Arenzon, A. J. Bray, L. F. Cugliandolo, and A. 

Sicilia, in preparation. 
[18] A. L. Stella and C. Vanderzande, Phys. Rev. Lett. 62, 

1067 (1989). W. Janke and A. M. J. Schakel, Phys. Rev. 

E 71, 036703 (2005). 
[19] D. Stauffer and A. Aharony, Introduction to percolation 

theory, 2nd ed. (Taylor & Francis, London, 1994). 
[20] A. J. Bray, K. Humayun, and T. J. Newman, Phys. Rev. 

B 43, 3699 (1991). 



